Atmos. Chem. Phys., 13, 9525-9541, 2013 
www.atmos-chem-phys.net/1 3/9525/20 1 3/ 
doi: 1 0.5 1 94/acp- 13-9525-2013 
© Author(s) 2013. CC Attribution 3.0 License. 



Atmospheric § 

Chemistry | MNCg 
and Physics | 


Aerosol airmass type mapping over the Urban Mexico City region 
from space-based multi-angle imaging 

F. Patadia 1 ’ 2 * , R. A. Kahn 2 , J. A. Limbacher 2,3 , S. P. Burton 4 , R. A. Ferrare 4 , C. A. Hostetler 4 , and J. W. Hair 4 

^ESTAR, Morgan State University, Baltimore, MD 21251, USA 

2 Climate & Radiation Laboratory, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 
3 Science Systems and Applications, Inc., Lanham, MD 20706, USA 
4 NASA Langley Research Center, Hampton VA, 23681, USA 

Correspondence to: F. Patadia (falguni.patadia@nasa.gov) 

Received: 6 January 2013 - Published in Atmos. Chem. Phys. Discuss.: 22 March 2013 
Revised: 30 July 2013 - Accepted: 31 July 2013 - Published: 27 September 2013 


Abstract. Using Multi-angle Imaging SpectroRadiome- 
ter (MISR) and sub-orbital measurements from the 2006 
INTEX-B/MILAGRO field campaign, in this study we ex- 
plore MISR’s ability to map different aerosol air mass types 
over the Mexico City metropolitan area. The aerosol air mass 
distinctions are based on shape, size and single scattering 
albedo retrievals from the MISR Research Aerosol Retrieval 
algorithm. In this region, the research algorithm identifies 
dust-dominated aerosol mixtures based on non-spherical par- 
ticle shape, whereas spherical biomass burning and urban 
pollution particles are distinguished by particle size. Two dis- 
tinct aerosol air mass types based on retrieved particle mi- 
crophysical properties, and four spatially distributed aerosol 
air masses, are identified in the MISR data on 6 March 
2006. The aerosol air mass type identification results are sup- 
ported by coincident, airborne high-spectral-resolution lidar 
(HSRL) measurements. Aerosol optical depth (AOD) gradi- 
ents are also consistent between the MISR and sub-orbital 
measurements, but particles having single-scattering albedo 
of ~ 0.7 at 558 nm must be included in the retrieval algorithm 
to produce good absolute AOD comparisons over pollution- 
dominated aerosol air masses. The MISR standard V22 AOD 
product, at 17.6 km resolution, captures the observed AOD 
gradients qualitatively, but retrievals at this coarse spatial 
scale and with limited spherical absorbing particle options 
underestimate AOD and do not retrieve particle properties 
adequately over this complex urban region. However, we 
demonstrate how AOD and aerosol type mapping can be 
accomplished with MISR data over complex urban regions, 
provided the retrieval is performed at sufficiently high spatial 


resolution, and with a rich enough set of aerosol components 
and mixtures. 


1 Introduction 

Aerosols play an important role in global climate forcing 
through their direct and indirect effects. They also have 

significant regional impacts on both climate and air qual- 
ity (IPCC, 2007). Information on aerosol type is crucial 
for improving our understanding and assessment of anthro- 
pogenic influences of aerosols on climate. Particulate mat- 
ter (PM), especially particles with aerodynamic diameter 
<2.5 pm, have known health impacts, such as cardiopul- 
monary diseases (Schwartz and Marcus, 1990; Saldiva et al., 
1995; Krewski et al., 2000). More generally, aerosol size in- 
formation is of value for particulate matter air quality assess- 
ment. The spatial and temporal distribution of aerosol mi- 
crophysical and optical properties is heterogeneous in na- 
ture and is caused by the wide variety of aerosol sources 
around the globe and relatively short atmospheric lifetime. 
In urban regions such as Mexico City, aerosols originate 
most commonly from industrial and domestic emissions, 
biomass burning, local transportation, and local wind-driven 
dust (Moffet et al., 2008; Stone et al., 2008; Aiken et al., 
2010). Aerosol sources and aerosol types can vary on short 

spatial scales (~ few km) in this region. Air quality stan- 
dards are frequently violated by high PM concentrations in 
the Mexico City area (Molina and Molina, 2004), and the 
aerosols also have the potential to affect the air quality and 
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climate of surrounding areas. For example, during the dry 
season, regional winds can transport aerosols from Mexico 
to the southeast United States (Wang et al., 2006; Voss et 
al., 2010). The study of aerosols therefore warrants contin- 
uous synoptic-scale monitoring, which can best be achieved 
through satellite measurements. 

Although satellites provide global monitoring of aerosol 
optical depth (AOD), few space-borne sensors have the ca- 
pability to distinguish different aerosol types. The Multi- 
angle Imaging SpectroRadiometer (MISR) aboard the Terra 
satellite is one instrument that can retrieve some informa- 
tion about aerosol shape, size and single scattering albedo 
(SSA) under good retrieval conditions (Kahn et al., 2005, 
2010). Particle type represents a constraint on particle ori- 
gin and composition, and can contribute to mapping aerosol 
evolution and transport. Under good retrieval conditions and 
depending on the aerosol type, the Version 22 MISR Stan- 
dard Aerosol Retrieval algorithm can distinguish two or three 
aerosol size groupings, about two SSA categories (absorbing 
vs. non-absorbing) and spherical vs. non-spherical particles. 
Over urban areas, Kahn et al. (2010) showed with MISR data 
that the aerosols are localized and vary over small spatial 
scales. This suggests the necessity for high-resolution AOD 
retrievals over urban areas. 

Comparisons of satellite aerosol property retrievals against 
suborbital measurements, such as those collected during field 
campaigns, provide unique opportunities to refine and assess 
the quality of the satellite retrievals. Using a combination of 
MISR standard aerosol retrievals (Martonchik et al., 1998, 
2009) and MISR research algorithm retrievals (Kahn et al., 
2001; Kahn and Limbacher 2012), along with ground and 
aircraft data from the Intercontinental Chemistry Transport 
Experiment-B (INTEX-B) field campaign, we analyze MISR 
aerosol retrievals in detail over the Mexico City area in this 
study. Although the MISR operational algorithm attempts to 
retrieve aerosol properties globally, and computational con- 
straints limit the range of aerosol models that can be consid- 
ered, the research algorithm allows us to explore the infor- 
mation content of the MISR data for a much broader range 
of particle types, and at significantly higher spatial resolution 
for this application. 

Our validation data comes from the host of instruments 
that obtained aerosol physical and optical property mea- 
surements during the INTEX-B field campaign. Our study 
adds information about the performance of MISR’s Standard 
aerosol retrieval algorithm to the existing literature, and also 
uses the satellite data to map AOD and aerosol type over 
the Mexico City Metropolitan area (MCMA) in more detail 
with the MISR research algorithm, providing regional con- 
text for the aerosol measurements made by multiple, subor- 
bital sensors. We focus on identifying different aerosol types 
in the complex Mexico City urban setting, and on mapping 
AOD gradients and aerosol air mass types to the extent pos- 
sible, using MISR data in the MCMA study region. Pre- 
vious studies relating to MISR aerosol-type sensitivity fo- 


cused on biomass burning and dust-dominated regions (Chen 
et al., 2008; Kalashnikova and Kahn, 2006; Kahn et al., 
2009a), and locations affected by thin cirrus (Pierce et al., 
2010), whereas this study focuses on the aerosols over the 
urban Mexico City basin. The paper is organized as follows: 
Sect. 2 briefly describes the datasets used, Sect. 3 discusses 
the sub-orbital constraints on aerosol properties for our study 
days, results from MISR research algorithm are discussed in 
Sect. 4, Sect. 5 presents notes on the sensitivity of MISR re- 
trievals over our study area and a summary of this study is 
given in Sect. 6. 


2 Study area and data used 

The Mexico City Metropolitan area, consisting of the heavily 
populated Mexico City basin, is our study region (see Fig. 1). 
Mexico City is one of the world’s largest megacities, and has 
undergone rapid urban development. We analyze data for two 
golden days, 6 and 15 March 2006, defined as days during 
which (1) MISR observations were made over the study re- 
gion, (2) coincident ground and air observations are avail- 
able, and (3) the key observational requirements of relatively 
cloud-free conditions and the presence of aerosols from dif- 
ferent sources are met. The ground-based and aircraft data 
used in this study are described briefly in this section, along 
with the Version 22 MISR Level 2 standard aerosol product. 

2.1 Ground-based and aircraft measurements 

During the field campaign, Aerosol Robotic Network 
(AERONET) CIMEL sun photometers were deployed at sites 
labeled as TO, T1 in Fig. 1, and a third at Tampico, about 
500 km northeast of Mexico City (not shown), and five Mi- 
crotops II sun photometers were deployed, one at each of the 
Hidalgo, UNAM, Corena, TEC, and UAMI sites (de Almeida 
Castanho et al., 2007; Molina et al., 2008), also shown in 
Fig. 1. The AERONET is a global, ground-based network of 
well-calibrated sun photometers that measure spectral AOD 
at 340, 380, 440, 500, 670, 870, 1020 and 1640 nm with an 
accuracy of about ±0.015 (Holben et al., 1998, 2001; Eck et 
al., 1999). AERONET also retrieves other column-effective 
aerosol properties, such as size distribution, refractive in- 
dex, single-scattering albedo (SSA), asymmetry parameter 
( g ) and phase function, with varying frequency and degrees 
of accuracy (Dubovik and King, 2000). We use version-2 
level-2 AOD and Angstrom Exponent (AE) data from the TO, 
T1 and Tampico stations for the two golden study days. The 
Microtops II are handheld sun photometers that measure so- 
lar radiance at five spectral wavelengths and provide AOD to 
within ±0.03 (Ichoku et al., 2002). We use the AOD at vis- 
ible wavelengths from all five ground-based Microtops sites 
to evaluate the corresponding values from both the MISR L2 
aerosol product and the MISR research algorithm. For each 
ground-based site, Table SI summarizes all the available 


Atmos. Chem. Phys., 13, 9525-9541, 2013 


www.atmos-chem-phys.net/13/9525/2013/ 


F. Patadia et al.: Aerosol airmass type mapping over the Urban Mexico City 


9527 



Fig. 1. (a) Nadir, 70° forward and 70° aft ward view of Mexico City 
from MISR (LIB RGB image) on 6 March 2006. Overlaid are loca- 
tions of ground-based instruments and the J3 1 (pink) and B200 (red) 
flight tracks (bottom right image). Names of the ground based sta- 
tions are indicated in the Figure. These also include the AERONET 
stations at supersites TO and Tl. The 70° forward view image shows 
fires burning at atleast 2 different locations on the mountains to the 
east of Mexico City and the transport of smoke from these fires over 
the MCMA study area, (b) Same as (a), except, it is for 1 5 March 
2006. More burning activity can be seen on this day (see top right 
image) compared to (a) for 6 March 2006. Also, notice the presence 
of clouds around TO supersite (top right image) and around central 
Mexico City (bottom right image). 


ground and aircraft instruments, the aerosol properties used, 
geographic location and the availability of different observa- 
tions on 6 and 15 March 2006. 

During the MILAGRO/INTEX-B field campaign, eight in- 
strumented research aircraft from the United States made 
measurements. Details on the principal airborne platforms 
that participated in the INTEX-B campaign and the scientific 
payloads on each aircraft can be found in Singh et al. (2009) 
and Molina et al. (2010). Rogers et al. (2009) describe the 
High Spectral Resolution Lidar (HSRL) carried by the B200 
aircraft. Only the J-31 and B200 aircraft acquired coinci- 
dent aerosol property measurements on both our case study 
days. For these days, we use primarily the vertical profiles of 
aerosol extinction, backscatter and depolarization, along with 


aerosol classification products, from the HSRL onboard the 
B200 to evaluate aerosol retrievals from the MISR standard 
and research algorithms. 

2.2 The MISR L2 aerosol product 

The Multi-angle Imaging SpectroRadiometer (MISR) is on 
board Terra, one of the NASA Earth Observing System 
(EOS) of satellites. Terra was launched on 18 December 
1999 into a sun-synchronous orbit at an altitude of 704 km, 
and has an equatorial crossing time of 10:30 a.m. LT (north 
to south). MISR images earth in four spectral bands centered 
at 446, 558, 672, 867 nm, and has nine push-broom cameras 
viewing at nine different angles, 0, ±26.1, ±45.6, ±60.0, and 
±70.5 degrees, covering the nadir, forward, and aft directions 
along the line-of-flight. The MISR level 2 aerosol data prod- 
uct (MIL2ASAE) contains: (1) aerosol optical depth (AOD) 
in 4 spectral channels, (2) fraction of AOD for three differ- 
ent particle size ranges, (3) fraction AOD assigned to non- 
spherical particles, (4) spectral single-scattering albedo, (5) 
angstrom exponent, other parameters related to aerosol re- 
trieval quality, and geographical information. In this paper, 
we use the MIL2ASAE Version 22 data product. The L2 
aerosol optical depth is reported at 17. 6x17. 6 km resolution. 
A detailed description of the MISR aerosol retrieval is given 
in Martonchik et al. (1998, 2009) and Kahn et al. (2009b). 

Validation of MISR AOD (558 nm) against AERONET 
shows that about 70-75 % of MISR AOD retrievals fall 
within 0.05 or 20 % of AERONET AOD and about 
50-55% are within 0.03 or 10% of AERONET AOD 
for globally distributed sites (Kahn et al., 2010). The 
AOD is designated as Stage 3 validated, whereas other 
aerosol properties such as SSA, AE and fractional 
AOD are Stage 2 validated, as described in the MISR 
Data Quality Statement distributed with the data prod- 
ucts (http://eosweb.larc.nasa.gov/PRODOCS/misr/Quality_ 
Summaries/misr qual stmts new.html). The MISR research 
algorithm used in this study is described in Sect. 3.3. 


3 Suborbital constraints on aerosol properties for the 
MISR golden days 

Figure la shows the MISR LIB RGB images (Nadir, 70° 
Forward and 70° Aftward cameras) from 6 March 2006 for 
our study region. The B200 flight track (with the HSRL in- 
strument) for the same day is superposed as a red line on the 
70° forward camera image. The J31 flight track is shown as 
a pink line on this figure, as this aircraft carried the AATS 
14 instrument (Redemann et al., 2009), from which a few 
measurements coincident with the MISR overpass were also 
obtained. The key ground-based sites are also marked on the 
same image. Figure lb is similar to Fig. la, except it is for 15 
March 2006. In the following section, we use the data from 
all these ground-based and airborne instruments to assess 
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aerosol characteristics, and spatial and temporal variability, 
on the two case study days. These results are then used to 
validate and test the limits of the MISR retrievals in the sub- 
sequent section. AOD from the J3 1 and the ground-based in- 
struments used in this study are interpolated to MISR mid- 
visible wavelength (558 nm) using the Angstrom formula, 
and is reported at 558 nm throughout this paper. 

3.1 Suborbital constraints on the aerosol environment 
for 6 March 2006 

On 6 March 2006, AOD is available over four of the Mi- 
crotops ground stations and over three AERONET stations 
(Tables Sib and Sic). In Fig. 2, the lines in different col- 
ors depict the AOD variations over different ground stations 
(mapped in Fig. 1) in the field campaign area. Aerosol opti- 
cal depth (Fig. 2a) shows diurnal variation with average AOD 
(558 nm) of 0.1 in the forenoon (16:00 UTC) and ~ 0.4 in the 
afternoon hours (by 21:00 UTC) over central Mexico City. 
Fower AOD observations spaced about 30 km apart at the 
T1 and Tampico sites to the northeast of Mexico City sug- 
gest gradients in AOD in the Mexico City basin. From the 
MISR RGB image for 6 March (Fig. la), we find at least two 
biomass-burning fires, one each to the northeast and south- 
east of Mexico City. The clockwise wind circulation (Fast 
et al., 2007) on this day explains the transport of aerosols 
from the fire locations into the field campaign region, pos- 
sibly missing the AERONET supersites farther north at T1 
and Tampico. Basin scale wind pattern on this day falls into 
the 03-South category (de Foy et al., 2008) that is character- 
ized by east-west convergence zones and southwest venting, 
which explains AOD pattern over the above AERONET su- 
persites. 

We can infer the type of aerosols in the study area by an- 
alyzing the vertical distribution of aerosol intensive proper- 
ties from the HSRF measurements (Burton et al., 2012). Fig- 
ure 3 shows the vertical distribution of aerosol properties on 
6 March 2006 as measured by the HSRF instrument on B200 
Flight 10. The HSRF made several observations around the 
Mexico City area from 16:75 to 17:75 UTC, and in the dis- 
cussion below we concentrate on this time window, as MISR 
passed over the area at 17:25 UTC. The aerosol extinction 
profiles at 532 nm (Fig. 3a) show high extinction at the sur- 
face and at 3.5 to 5 kma.s.l. (see red colors). This shows the 
high aerosol concentrations both in the boundary layer and 
higher in the atmosphere. During the MIFAGRO campaign, 
the maximum daily mixing layer depths always reached 

4.2 km a.s.l. and frequently extended to 6.2 km a.s.l. (Shaw et 
al., 2007). Mexico City is at an elevation of ~ 2.24 km a.s.l. 
Figure 3b, c and d, respectively, present the aerosol extinc- 
tion to backscatter ratio (532 nm), aerosol wavelength depen- 
dence (1064 nm/532 nm) or the backscatter related Angstrom 
exponent, and the aerosol depolarization ratio (532 nm). The 
aerosol extinction to backscatter ratio (Fig. 3b) has high val- 
ues near the surface at some locations (e.g., see values corre- 


06-1V larch-2006 




Fig. 2. (a) Diurnal variation of AOD (558 nm) on 6 March 
2006, over the ground based stations (shown in Fig. 1) in 
MILAGRO/INTEX-B field campaign area around Mexico City. 
Aerosol optical depth from AERONET stations at supersites TO, 
T1 and Tampico are also shown, (b) Same as (a) but for 15 March 
2006. 


sponding to 17:20-17:25 UTC). In this time window, aerosol 
wavelength dependence (1064/532 nm) shows a similar pat- 
tern (Fig. 3c) whereas the aerosol depolarization ratio is low 
(Fig. 3d). High extinction-to-backscatter ratio, high wave- 
length dependence and low depolarization values can indi- 
cate the presence of urban aerosol, as shown in the HSRF 
classification result in Fig. 3e. On the other hand, high de- 
polarization values (e.g., below ~ 3 km and between 17:35— 
17:40 UTC in Fig. 3d) are indicative of the presence of dust 
particles (see Fig. 3e) (Burton et al., 2012). Qualitative classi- 
fication of aerosol intensive parameters measured by HSRF 
(Burton et al., 2012), shown in Fig. 3e, indicates the pres- 
ence of three distinct aerosol types (dusty mix, smoke and ur- 
ban pollution particles) on 6 March. It is noteworthy that the 
dusty mix type identified by HSRF here is not pure dust. This 
identification occurs when depolarization measurements in- 
dicate the presence of dust along with other components (par- 
ticularly smoke or pollution) that are not distinguishable by 
the HSRF (see Burton et al., 2012 for details). Several other 
studies that analyzed the aerosol optical and chemical proper- 
ties during the field campaign (e.g., Marr et al., 2006; Marley 
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Fig. 3. HSRL observation on 6 March 2006: (a) aerosol extinction (b) aerosol extinction to back-scatter (Lidar) ratio (c) aerosol wavelength 
dependence (d) aerosol depolarization ratio (e) vertical distribution of aerosol types identified by HSRL (f) AOD attributed to aerosol types 
identified by HSRL. The approximate intersections of HSRL observations with the four aerosol air masses, as identified by MISR, are marked 
in (e) and (f). The legend for the air mass is shown in (f). Note that air masses over some region were observed multiple times by HSRL as 
the B200 aircraft circled in the study region (see Fig. la for B200 flight path). The elevations shown are above sea level (a.s.l.). 


et al., 2007, 2009a; Stone et al., 2008; Querol et al., 2008; 
Corr et al., 2009; Paredes-Miranda et al., 2009; Bergstrom et 
al., 2010; De Foy et al., 2011) corroborate the aerosol type 
results from the HSRL and add detail, suggesting that the 
study area contained urban aerosols, including black carbon 
and organic matter, and mineral dust aerosols during the cam- 
paign period. Ground-based instruments also observed three- 
to-four aerosol size modes that are likely from these different 
aerosol types in the study area. 

3.2 Suborbital constraints on the aerosol environment 
for 15 March 2006 

In contrast to 6 March, the AOD from ground stations is 
higher in the study area on 1 5 March and there is not much 
variation in the AOD throughout the day (Fig. 2b). How- 
ever, similar to 6 March, a spatial AOD gradient between 
central Mexico City and its surrounding region is observed 
in Fig. 2b. The lack of diurnal variation on 15 March can 


be attributed to locally uniform spatial mixing of aerosols 
due to high winds (15 ms -1 ) associated with a cold surge 
event in the study area (Fast et al., 2007). Cold Surge events 
vent in the southwest to north to southeast direction, with 
less southward venting (de Foy et al., 2008). Kleinman et 
al. (2009) analyzed the DMA and PCASP aircraft measure- 
ments of aerosol size distribution and found three-to-four 
modes in MCMA on 15 March. This aerosol size distribu- 
tion is produced by multiple aerosol sources in the region. 
Moffet et al. (2008) and Stone et al. (2008) suggest that the 
submicron-mode aerosols are comprised of industrial emis- 
sions, biomass and biofuel burning and aged organic carbon 
(OC). Aiken et al. (2010) show that the dominant source of 
BB OC aerosols found in the Mexico City area are from for- 
est fires burning in the surrounding mountains. The coarse 
mode was dominated by inorganic dust particles, 50% of 
which were coated with OC (Moffet et al., 2008). Wang et 
al. (2010) suggest traffic emissions, aging of aerosols, lo- 
cal photochemical production of secondary aerosol species 
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Fig. 4. HSRL observation on 15 March 2006: (a) aerosol extinction (b) aerosol extinction to back-scatter (Lidar) ratio (c) aerosol wavelength 
dependence (d) aerosol depolarization ratio (e) vertical distribution of aerosol types identified by HSRL (f) AOD attributed to aerosol types 
identified by HSRL. Blue box in (e) and (f) depicts the flight path over Mexico City. 


and wind-blown dust as other sources of aerosols in MCMA. 
Analysis of the HSRL-measured intensive aerosol properties 
over our study area between 17:05 and 17:35 UTC on this 
day (Fig. 4) also shows that urban aerosols accounted for 
nearly 80 % (median value) of the total AOD. 

Based on the analysis of aerosol properties measured by 
ground-based and aircraft instruments on our case study 
days, the aerosol types found in the study area are again 
comprised of smoke, urban pollution and dust aerosols. How- 
ever, the meteorology of the region is different on these days, 
which leads to differences in the spatial and temporal distri- 
butions of aerosols. Wang et al. (2010) suggest that during 
the daytime, over the Mexico City area, changes in aerosol 
size and composition occur over time scales of about a few 
hours. Spatial gradients in AOD and size distribution mea- 
sured at various surface sites also suggest aerosol variability 
at small spatial scales. The selection of aerosol models for 
MISR research retrievals (discussed in Sect. 4) is based on 
published literature and the analysis of suborbital informa- 
tion on aerosol type presented in Sects. 3.1 and 3.2. 


4 The MISR aerosol research retrievals with inputs 
from field campaign data 

4.1 Research retrieval details 

4.1.1 Research retrieval theory 

In the MISR Research Aerosol Retrieval algorithm (Kahn et 
al., 2001; Kahn and Limbacher, 2012), the TOA radiation is 
simulated using a radiative transfer (RT) code developed by 
the MISR Team (Martonchik et al., 1998) that is based on 
the matrix operator method (Grant and Hunt, 1968). For a 
given aerosol particle analog, vertical distribution, sun and 
viewing geometry and surface reflectance properties, the RT 
code simulates the TOA radiances as would be observed by 
the MISR instrument. The radiances produced by up to four 
component aerosol particles are then weighted by their frac- 
tional contributions to the AOD to obtain the radiance of 
the aerosol mixture (Abdou et al., 1997). The over-land re- 
trievals use all four MISR wavelengths. In the present study, 
MISR observations are simulated over MCMA. We assume 
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a Rayleigh-scattering atmosphere with a surface pressure of 
776 mb (Mexico City is at an elevation of ~ 2200 m), and 
a standard mid-latitude temperature profile. The ozone abun- 
dance in the study area was obtained from Ozone Monitoring 
Instrument (OMI) data (Ahmad et ah, 2003). The aerosols are 
assumed to be concentrated in a near-surface layer, except for 
medium-mode dust aerosols. The surface reflectance param- 
eters are prescribed from the MISR Land Surface Product at 
each of the four MISR wavelengths, based on an analysis of 
the MISR surface retrieval record on low-AOD days, for the 
same viewing geometry and location of interest. In this case, 
the atmosphere is relatively cloud-free, the surface is reason- 
ably well characterized, and the AOD is sufficient to provide 
a good atmospheric signal. These are reasonably good atmo- 
spheric conditions for remote-sensing aerosol retrievals. 

The TOA simulated radiances are converted to TOA re- 
flectance and compared to the MISR observed reflectance 
to determine the AOD and aerosol type. The agreement be- 
tween simulated and observed TOA reflectance at all four 
wavelengths and at 9 viewing angles is assessed using three 
different Chi-Squared (x abs , Xg eom > X S p CC >) tests ( Kahn et al -> 
2001). The largest of (x abs , Xg e0 m> X^pec) for a § iven mixture 
is designated as X,, lax _3 - The comparisons are made for all 
possible combinations of aerosol particles in different pro- 
portions (from 0 to 100 % in increments of 5 %) within each 
four-component mixing group. In our study, the AOD for 
each combination is varied from 0.0 to 1.0, in steps of 0.05. 
We examined a hierarchy of relative and absolute acceptance 
criteria applied to ChiMin (the minimum value of Xmax-3 
for all mixtures in a retrieval region), and looked for con- 
vergence toward a consistent result as more stringent crite- 
ria were applied. Kahn and Limbacher (2012) used similar 
criteria and provide further details on this approach. Unless 
specified otherwise, values shown in this paper are based 
on those retrievals meeting the Xmax-3 < ChiMin + 0.1 cri- 
terion. For a given HR region, the mean values of aerosol 
properties (e.g., total AOD, small/medium/large mode and 
spherical/non-spherical AOD fractions, SSA, AE) from all 
mixtures that satisfy our Xmax-3 acce P tance criteria are re- 
ported as successful retrievals and used in our analysis. 

4.1.2 Aerosol components and mixing groups 

Table 1 presents the four-component aerosol mixing groups. 
Table S2 lists the aerosol component particles that define 
the aerosol mixtures used in this study. Each of the 45 dis- 
tinct aerosol components in Table S2 is defined by spec- 
tral real and imaginary indices of refraction and a lognormal 
size distribution, given as the particle geometric mean radius 
( r pg,N) and geometric standard deviation (a g ). The optical 
properties for each spherical component are derived using a 
standard Mie scattering code, with the exception of the dust 
grains. The component particles represent biomass burning, 
urban pollution and dust aerosol particles that span five par- 
ticles size ranges (very very small (VVS), very small (VS), 


small (S), medium (M) and large (L) particles) and four sin- 
gle scattering albedo ranges (very absorbing (VA), absorb- 
ing (A), less absorbing (LA) and weakly absorbing (WA) 
particles). For biomass burning particle analogs in Table S2, 
the sizes are: VVS (r pg? n = 0.04 pm), VS (r pg? n = 0.07 pm), 
S (r pg ,N = 0.13 pm) and S2 (r pg ,N = 0.16 pm) and the mid- 
visible SSA values are: A (<u 0 558 = 0.84), LA (&> 0 558 = 0.90) 
and WA (<u 0 558 = 0.94). The sizes of the other spherical par- 
ticle types that represent both pollution and biomass burning 
particle analogs (Table S2) are: VVS (r pg? n = 0.06 pm), VS 
(r pg N = 0.12 pm), S (r pg? N = 0.24 pm), M (r pg? N = 0.24 pm) 
and L (r pg? n = 0.5 pm) and the mid- visible SSA values 
are: VA (&> 0 558 = 0.8), A (<u 0 558 = 0.85), LA (<u 0 558 = 0.9) 
and WA (<u 0 558 = 0.95). Strong spectral dependence of 
aerosol absorption defines the “steep” particles that repre- 
sent biomass burning aerosols, whereas the “flat” particles 
that represent pollution aerosols have spectrally independent 
absorption (Bond and Bergstrom, 2006). 

Table 1 lists the 23 mixing groups used in this study. These 
mixtures are created to test the ability of the MISR algorithm 
to distinguish between sizes (M1-M4), steep and flat SSA 
spectral slope (M6, M7, M21), spherical and non-spherical 
(M10-M13), absorbing and non-absorbing (M14-M17) par- 
ticles (Kahn et al., 2001). Based on the targeted tests, the 
mixing groups are collected into 10 larger groupings listed in 
column 1 of Table 1 . The last column in Table 1 provides a 
detailed description of the objective for each grouping. The 
mixing groups are created such that they span the range of 
sizes and optical properties of biomass burning, pollution 
and dust particles that were found in the MCMA during the 
INTEX-B field campaign. 

4.1.3 Research retrieval regions 

Figure 5 shows the area over which MISR research re- 
trievals are performed in this study for 6 March. There 
are more than 400 individual retrieval regions in the re- 
trieval area in Fig. 5. Research retrievals are performed over 
825 m x 825 m areas (3 x 3 275 m MISR Local Mode radi- 
ance pixels) and this defines the retrieval regions in our case. 
We call these “high-resolution” retrievals when compared to 
the 17.6 km x 1 7.6 km retrieval regions of the L2 MISR Stan- 
dard aerosol product, and refer to them as HR regions in the 
paper. Note that the MISR instrument can be configured to 
provide high-resolution (275 x 275 m) images in all 36 chan- 
nels for targets and observation times selected several weeks 
in advance of the overpass. This capability is referred to as 
Local Mode. In the default Global Mode, all nadir cameras 
and all off nadir cameras for Red channels (Band 3) are kept 
at high resolution (275 m) whereas observations in the re- 
maining 24 cameras are down-sampled onboard to 1.1 km 
(http://www-misr.jpl.nasa.gov/getData/localMode/). The lo- 
cations of ground-based Microtops and AERONET stations 
are indicated as small white circles in Fig. 5. If ground ob- 
servations of AOD are available at satellite overpass time, 
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Table 1 . Aerosol Mixing Groups. 


Groups 

Mixing 

Groups 

Component 1* 

Component2 

Component3 

Component4 

Mixture Description 


Ml 

pl9 

p29 

p33 

p37 

These 4 mixtures have particles with a constant SSA but 

Group 1 

M2 

pl6 

p20 

p30 

p34 

different particles sizes. This group of mixtures (Ml- 

M3 

P2 

p4 

P7 

plO 

M4) tests the sensitivity of MISR algorithm to particle 


M4 

p38 

p39 

p40 

p41 

size for our study and the most likely SSA in the study 
area. The sizes are VVS, VS, S, M in the 4 mixtures 
and the SSA of each mixture is LA, WA, LA, and Non- 
absorbing respectively. 

Group 2 

M5 

p31 

p32 

p35 

p36 

This mixture of VA and A type of S and M steep parti- 
cles tests if a combination of size and SSA yields better 
results compared to above group of mixtures with con- 
stant SSA 

Group 3 

M6 

pi 7 

pl8 

pl9 

p20 

Mixture of Flat Small or Medium particles with differ- 


M7 

p21 

p22 

p23 

p24 

ent SSA is used to test the sensitivity to distinguish pol- 
lution from smoke particles 

Group 4 

M8 

pl3 

p27 

p29 

P 15 

These mixing groups test the likely SSA of parti- 


M9 

pl4 

p20 

p29 

p33 

cle sizes that produce acceptable results. Sensitivity to 
steep and flat particles is also tested. 

Group 5 

M10 

pl8 

p22 

p32 

P 43 

Same as the above group of mixtures but with a dust 


Mil 

pl9 

p33 

p45 

— 

component (2 different sizes) added to the flat and steep 
particles of different sizes (S, M). 

Group 6 

M12 

pl9 

p33 

p37 

p43 

This group tests if dust can be replaced by a non- 


M13 

p25 

p29 

p33 

p41 

absorbing medium sized particle. Ml 2 consists of dust 
and Ml 3 has a medium non-absorbing particle. 

Group 7 

M14 

pl2 

pl5 

p33 

p41 

This group of mixtures consists of particles with differ- 


M15 

p25 

p29 

p42 

p43 

ent sizes and SSA and tests MISR algorithm’s sensitiv- 


M16 

pl9 

p40 

p41 

p42 

ity to absorbing vs non-absorbing particles. 


M17 

P 33 

p37 

p42 

p43 


Group 8 

M18 

p2 

p4 

P7 

p9 

This group of mixtures consists of biomass burning 


M19 

p2 

p4 

p6 

plO 

(BB) particles used by Chen et al. (2008). They are 


M20 

p3 

p5 

P 8 

p26 

used to test if the retrieval algorithm favors these par- 
ticles over their spherical absorbing counterparts used 
in groups 1-7. The sensitivity to size and SSA of BB is 
also tested in this group. 

Group 9 

M21 

p2 

p6 

pll 

pl8 

Mixture of BB and pollution like VVS and S particles 
with LA and A SSA. Similar to group 3, this group tests 
the algorithm’s sensitivity to distinguish pollution from 
smoke particles 

Group 10 

M22 

P2 

p9 

p44 

- 

Same as the group 5 but BB particles. 


M2 3 

Pi 

p2 

P? 

p43 



* The aerosol particle components (pl-p45) in columns 3-6 are described in Table S2. 


their locations are indicated as red circles in Figs. 5a and 6a 
for 6 and 15 March, respectively. The HR region closest to 
the ground station having a valid AOD observation is indi- 
cated by a red-black circled dot in Figs. 5a and 6a. HSRL 
AOD along the flight track is also shown for 6 and 1 5 March 
in Figs. 5a and 6a, respectively. AOD comparisons are tabu- 
lated in Table 2. The following section discusses the MISR 
research retrieval results. 


4.2 MISR research retrievals over Mexico City Basin on 
6 March 2006 

An analysis of the AOD and aerosol properties (shape, size, 
and SSA) from the MISR research retrieval algorithm, iden- 
tifying aerosol air mass type, is presented here. Unless spec- 
ified otherwise, the aerosol air mass and associated prop- 
erties refer to MISR Research Algorithm retrievals. Fig- 
ure 5a shows the spatial distribution of 558 nm AOD from the 
MISR research algorithm, overlaid with 558 nm AOD from 
ground-based observations and 532 nm HSRL AOD. AOD 
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Fig. 5. Aerosol properties from MISR research retrieval algorithm over MCMA for 6 March 2006: (a) aerosol optical depth at 558 nm The 
HSRL 532 nm AOD along the flight path is superposed in panel (a). For HSRL, black color indicates missing data, (b) Non-spherical fraction 
of total AOD (c) Single Scattering Albedo at 558 nm (d) fraction of AOD by 0.06 pm particles (e) fraction of AOD by 0.12 pm particles (f) 
fraction of AOD by 0.2-0.28 pm particles, (d, e, f) represent the small mode fraction of AOD. Distribution of large mode fraction of AOD is 
similar to (b). The 4 aerosol air masses identified by MISR are shown in (c, d, e, f). 


from available coincident measurements are tabulated in Ta- 
ble 2. MISR AOD agrees to within 0.03 of Microtops AOD 
at UAMI and TEC. However, MISR AOD (0.13) is about 
half that at the TO supersite (0.25) in the central city. Given 
the low AOD (0.25) condition on this day, MISR research 
retrievals show ambiguity in the choice of particle size and 
SSA over the TO site. This results in a large number of ac- 
ceptable AOD retrievals that range from 0.05-0.3, the mean 
of which is 0.13. The higher AOD retrievals, however, cor- 
respond to mixtures comprising of VVS, VS and S particles 
and lower SSA (< 0.80). HSRL AODs (532 nm) are within 
0.05 of all ground-based AODs (532 nm shown in bold in 
Table 2). Note that the HSRL extinction measurement starts 
at 300 m above ground level. The extinction from 300 m to 
ground is extrapolated assuming constant mixing ratio in the 
boundary layer, and then used to calculate the total column 
AOD. In this paper (e.g., in Table 2) we call this the HSRL 
AOD. AOD is also calculated at each of the vertical layers 
(30 m) using HSRL extinction measurements. We call the in- 
tegrated layer AOD (without the AOD extrapolated below 
300 m) as LAOD in this paper. The LAOD is apportioned 
into different aerosol classes. 

Ligure 5b shows the spatial distribution of the fraction of 
non-spherical particles in the aerosol mixtures that produce 
successful MISR AOD retrievals. The spatial distribution of 


SSA and fraction of AOD attributable to 0.06, 0.12 and 0.22- 
0.26 pm sized particles in the aerosol mixtures are shown 
in Lig. 5c, d, e and f, respectively. The regions marked as 
I and II in Lig. 5b are dominated by spherical absorbing par- 
ticles (analogs of smoke and pollution particles) uniformly, 
whereas other regions in the retrieval area consist of both 
spherical and non-spherical (dust analog) particles. The other 
aerosol properties (AOD, SSA and sizes) are less uniform in 
regions I and II and elsewhere. We identify four spatially 
distributed aerosol air masses that appear to have distinct 
aerosol properties in the MISR retrievals and for which coin- 
cident HSRL observations are available for evaluation. These 
aerosol air masses are marked in Lig. 5c-f, and their char- 
acteristics from MISR and HSRL are tabulated in Tables 3 
and 4, respectively. According to MISR, aerosol air mass Al, 
which is downwind of the NE fires, consists predominantly 
of spherical, absorbing smoke particles from biomass burn- 
ing. Lrom Table 3 we see that A2, A3 and A4 are similar in 
terms of their SSA range and the percent contribution from 
0.12 and 0.26 size particles, and the differences are within 
the range of measurement uncertainties. Aerosol type clas- 
sifications from HSRL (see Table 4) also show more smoke 
(~ 55 %) in Al and more urban pollution particles in A2, 
A3 and A4 (41-52%). HSRL data imply that the prevalent 
0.26 pm particles observed by MISR in Al are analogs of 
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Table 2. AOD (at 558 nm and 532 nm) comparisons between MISR Level 2 (both standard and research retrievals) and other available ground 


based and aircraft instruments. 

Date Site 

Instrument 

Time a 

(UTC) 

Distance 
from Site b 

AATS-14 

Pressure 

Altitude 

Sub-orbital 
558 AOD (532 
AOD) 

MISR-R 

MISR L2 

6- Mar-0 6 

UAMI 

Microtops 

17.33 

- 

- 

0.21 (0.22) 

0.18 (778 m) 



HSRL C 

17.10 

2.49 km 

- 

-(0.19) 




AATS-14 

17.28 

2.74 km 

2.63 km 

0.13(0.15) 




CIMEL 

17.32 

- 

- 

0.25 (0.27) 



TO-AERONET 

HSRL 

16.84 

1.41km 

- 

-(0.21) 

0.13 (550 m) 

0.10 (9 km) 


AATS-14 

17.47 

230 m 

2.63 km 

0.11 (0.13) 




Microtops 

17.33 

- 

- 

0.13 (0.15) 



TEC 

HSRL 

- 

- 

- 

- 

0.11 (615 m) 

0.1 (2.9 km) 


AATS-14 

16.8 

5.65 km 

5.35 km 

0.005 (0.005) 



15-Mar-06 

UAMI 

Microtops 

17.25 

- 

- 

0.40 (0.43) 

0.41 (2.2 km) 

0.22 (6.7 Ion) 


HSRL 

17.18 

4.2 km 

- 

- (0.44) 




AATS-14 

17.86 

6km 

4.9 km 

0.25 (0.25) 




Microtops 

17.25 

- 

- 

0.54 (0.58) 



HIDALGO 

HSRL 

17.16 

1.63 

- 

(0.47) 

0.28 (1.13 km) 

0.12 (5.4 km) 


AATS-14 

17.25 

922 m 

2.6 km 

0.19(0.20) 



UNAM 

Microtops 

- 

- 

- 

- 

0.43 
(2.3 Ion) 

- 


HSRL 

17.22 

2.3 km 

- 

- (0.32) 



CORENA 

Micro tops 

17.7 



2.5 

(anamolous) 

- 

- 


HSRL 

17.36 

1.75 km 

- 

0.32 




CIMEL 

17.48 

- 

- 

0.41 (0.44) 



TO-AERONET 

HSRL 

17.14 

42.6 m 

- 

-(0.43) 

0.25 (1.85 Ion) 

0.12 (9.4 Ion) 


AATS-14 

17.23 

112m 

2.63 km 

0.22 (0.23) 



a MISR overpass time is 17:37 UTC on 6 March 2006 and 17:26 UTC on 

15 March 2006. 






b Distance of the aircraft measurement from the surface Site listed in column 2. In the last two columns, values in parentheses are distances from the surface site to the 
center of the spacecraft retrieval region. 

c The HSRL provides layer AOD (~ 0-~ 7 km AGL) unlike columnar measurements from sun photometers. HSRL AOD is at 532 nm. 

No AATS-14 observation over TEC ground site. 

MISR-R is the AOD from MISR research retrieval algorithm averaged over all mixtures that pass the acceptance criteria, and correspondingly, MISR L2 is the standard 
MISR “best estimate” AOD product. 

Note: sub-orbital AODs at 558 nm in column 7 are calculated using the angstrom formula. 


smoke particles and the smaller 0.12 pm particles in A2 are 
better analogs of pollution particles. MISR retrievals over the 
central Mexico City area show that the aerosol air mass A3 
is dominated by pollution particles (0.12 pm) that are mixed 
with ~ 35 % dust. HSRL observations verify the presence of 
more urban (51 %) aerosols in A3 compared to Al, although 
the percent contribution of dusty mix to LAOD is similar 
(Table 4). Ground-based instruments that observed inorganic 
dust particles and sub-micron sized aerosols in the A3 area 
provide further validation of dust detected by MISR in A3. 
The broader size distribution of small particles in MISR re- 


trievals over Al, A2, and A4, however, is not supported by 
the HSRL data, which instead shows ~ 39-45 % contribution 
from dusty mix type aerosols. However, it is noteworthy that 
HSRL may have only observed a portion of the aerosol air 
masses identified by MISR e.g., HSRL observed the right- 
most edge of A3 where MISR does see up to 45 % dust. Also 
note that the dusty mix HSRL classification is not pure dust 
and this identification occurs when depolarization measure- 
ments indicate the presence of dust with other components 
(particularly pollution vs. smoke) (Burton et al., 2012). Ta- 
ble 4 shows the uncertainty in AOD assigned to different 
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Table 3. Characteristics of aerosol air masses based on shape, size and single scattering albedo as identified from the MISR research retrievals 
over the Mexico City metropolitan area. 



Air Mass 
(size) a 

Shape 

AOD 
(558 nm) 


Size (r e ff 

in pm) %Contribution 


SSA 

(558 nm) 


Date 



Range 

0.06 

0.12 

0.22-0.26 

0.57 

0.75 


Comment 

6- Mar-0 6 

Al 

(~ 75 km 2 ) 

Spherical 

0.25-0.35 

<25 

10-30 

30-80 

- 

- 

0.85- 

0.9 

Region comprised of small 
mode particles (reff < 0.26 pm) 


A2 

(~ 50 km 2 ) 

Spherical 

0.1-0.25 

5-35 

15-50 

35-45 

<20 

— 

0.9- 

0.95 

Region comprised of small (reff 
< 0.26 pm) and medium (reff = 
0.57 pm) mode particles 


A3/MC b 

(~35km 2 ) 

Spherical 
+ Non- 
spherical/ 
Dust 

0.1-0.18 

<25 

15-50 

25-35 


10-40 

0.9- 

0.95 

Region comprised of small (reff 
< 0.26 pm) and large (reff = 
0.75 pm) mode particles 


A4 

(~25 km 2 ) 

Spherical 

0.1-0.18 

<25 

20-40 

15-35 

<30 


0.9- 

0.95 

Region comprised of small (reff 
< 0.26 pm) and medium (reff = 
0.57 pm) mode particles 

15-Mar-06 

NE Quad 

Spherical 

0.2-0.45 

<25 

10-50 

10-45 

5-55 


0.9-1. 0 

Region comprised of small (reff 
<0.26 pm) and medium (reff = 
0.57 pm) mode particles 


NW Quad 

Spherical 

0.2-0.35 

5-40 

5-50 

10-45 

<40 

— 

0.9- 

0.95 

Region comprised of small (reff 
< 0.26 pm) and medium (reff = 
0.57 pm) mode particles 


SE Quad 

Spherical 

0. 1-0.3 

15-60 

<25 

<30 

<50 


0.88- 

0.95 

Region comprised of small 
(reff < 0.26 pm) and (reff 
<0.26 pm)mode particles 
Note that the AODs closer to 
Mexico City are cloud contam- 
inated 


SW Quad 

Spherical 
+ Non- 
spherical/ 
Cirrus 

Cloud Con- 
taminated 







Cirrus clouds found in this area 


a This is the approximate size of the aerosol air mass. The aerosol air mass sizes were estimated based on visual reading of the number of pixels (that fell within each air 
mass we identified) multiplied by size of single retrieval pixel. 
b MC = Mexico City. 

Note: numbers in bold highlight the particle sizes in each mixture that makes highest % contribution to the total AOD. 


Table 4. Characteristics of aerosol air masses from HSRL data over the Mexico City metropolitan area. These air masses correspond to the 
air masses identified by MISR in Table 5. 


Air mass 

HSRL- AOD Range 
over each air mass 

%DUST- AOD 
(Uncertainty) 

% Smoke- AOD 
(Uncertainty) 

%Urban-AOD 

(Uncertainty) 

Al 

0.14-0.30 

37% (7%) 

55% (10%) 

4 % (1 %) 

A2 

0.21-0.25 

45% (19%) 

3% (1.5%) 

52% (22%) 

A3 

0.19-0.25 

38% (15%) 

10% (4%) 

50% (20%) 

A4 

0.18-0.25 

39% (9%) 

14% (3%) 

41 % (9 %) 


Note: (1) The DUST in column 3 pertains to both “dusty mix” and pure dust. The dusty mix identification occurs 
when depolarization measurements indicate the presence of dust with other components (particularly pollution vs. 
smoke) (Burton et al., 2012). (2) Smoke pertains to both smoke and fresh smoke classes (3) % Contribution of 
Smoke, Dust and Urban AOD values are weighted mean values (weighted by AOD). The percentage in brackets 
indicate the uncertainty in %contribution for the given aerosol type. This is the uncertainty due to the fraction of 
total AOD that remains unclassified from extinction measurements below 300 m. 


aerosol classes. We calculate the uncertainty value by consid- 
ering HSRL AOD as the total AOD instead of the LAOD; this 
accounts for unclassified near-surface AOD and provides the 
lower limit on the calculations. However, as this location is 
a source region for pollution particles, the near-surface com- 
ponent is likely to contain a larger fraction of such particles 


than the total column. Also, note that MISR is generally not 
sensitive to the properties of aerosol components that con- 
tribute < 20 % to the total AOD (Kahn et al., 2001), which in 
this case could include the dust component reported by the 
HSRL over A2 and A4, where MISR detects up to 20 and 
30 % medium sized particles but no dust. 
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Fig. 6. Aerosol properties from MISR research retrieval algorithm over MCMA on 15 March 2006: (a) aerosol optical depth at 558 nm. The 
HSRL 532 nm AOD along the flight path is superposed in panel (a). For HSRL, black color indicates missing data (b) non-spherical fraction 
of AOD (c) Single Scattering Albedo at 558 nm (d) fraction of AOD by size 0.06 pm (e) fraction of AOD by size 0.12 pm (f) fraction of AOD 
by size 0.26 pm. (d, e, f) represent the small mode fraction of AOD. Distribution of large mode fraction of AOD is similar to (b). 


The comparison of aerosol properties from MISR against 
sub-orbital HSRL observations illustrates the capability of 
MISR to qualitatively distinguish aerosol air masses in the 
MCMA, and suggests that the relative differences in small 
vs. coarse mode particles in the four aerosol air masses are 
robust. From the comparisons with HSRL we find that on this 
day, of the four air masses (A1-A4), MISR is able to distin- 
guish two air mass types in the study region based solely on 
particle microphysical properties - one dominated by smoke 
particles (Al) and the others that are mixtures of smoke, pol- 
lution and/or dust particles (A2-A4). And as expected, the 
aerosol type information from MISR, although significant, 
must be viewed as qualitative. 

4.3 MISR research retrievals over Mexico City basin on 
15 March 2006 

Figure 6a shows the mid- visible AOD retrieved by the MISR 
research algorithm in our study region on 15 March 2006. 
AOD ranges between 0.15 and 0.5, higher than that on 6 
March (Fig. 5a). The MISR LIB true-color image in Fig. lb 
shows more smoke from fires surrounding the Mexico City 
basin on this day than on 6 March (compare Fig. la, b), 
and a mixture of clouds and fires in the western half of the 
study region. Note that MISR cannot perform aerosol re- 
trievals over cloud-covered regions. A gradient can be seen 


in both AOD and MISR particle property maps in Fig. 6. In 
the northern half of the figure, there is an east-west gra- 
dient in the AOD that is likely due to varying concentra- 
tions of biomass burning aerosols from the fires on the NE 
side of the Mexico City basin. Northerly winds as high as 
15 ms -1 were observed in the area on 15 March (Fast et al., 
2007). Gradients in AOD in the SE quadrant also corrobo- 
rate the gradient in smoke concentration seen qualitatively 
in Fig. lb. Where available, the retrieved AOD comparison 
to ground-based AOD is shown in Fig. 6a (red circles). At 
MISR overpass time (17:26 UTC), a gradient in AOD is also 
reflected in the surface measurements, with higher AOD at 
Hidalgo (0.54) than at the TO (0.4) or UAMI (0.41) sites. Hi- 
dalgo is closer to Mexico City, whereas TO is to the north, 
and UAMI is to the southeast, as shown in Fig. 6, providing 
qualitative verification of the MISR-retrieved AOD pattern. 
Near-coincident HSRL observations provide further valida- 
tion, with higher AOD (0.47) in the city (near Hidalgo) and 
lower AOD outside the city (AOD = 0.43 at TO and 0.32 at 
TEC and UNAM/Corena) (see Table 2 for AOD compar- 
isons). At UAMI, the MISR AOD (0.41) is in good agree- 
ment with the Microtops and HSRL AODs (0.40, 0.44 re- 
spectively), whereas at other ground sites, MISR AOD is 
cloud contaminated (see Fig. lb. Notice the differences in 
the locations of clouds relative to the surface in the forward 
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and aft camera views that resulted in the cloud-contaminated 
retrievals). 

Returning to Fig. 6 b, the fraction AOD non-spherical is 
higher (50-80%) around central Mexico City and lower 
(generally < 10 %) elsewhere. This is due to the presence of 
the thin patchy cirrus clouds at high altitudes (> 9.5 km a.s.l.) 
in the region, as confirmed by the J3 1 and B200 aircraft field 
catalogs (http://catalog.eol.ucar.edu/cgi-bin/milagro/report/ 
index). Examination of cirrus reflectance from MODIS also 
indicates the presence of thin localized cirrus clouds over the 
central city. However, the HSRL aerosol classification in the 
SW quadrant indicates that urban aerosols made the great- 
est contribution to AOD over the city, whereas a mixture of 
dust and urban aerosols were the dominant retrieved AOD 
components outside the city (Fig. 4). In other parts of the 
study area (NW, NE and SE quadrants), the particle proper- 
ties agree qualitatively with the presence of smoke, as ob- 
served by MISR (see Fig. lb). The research algorithm favors 
~ 50 % VVS and ~ 30 % VS sized biomass burning aerosols 
in the NE quadrant. The HSRL transect to the north also 
shows >55 % smoke in this region (classification in this re- 
gion not shown in Fig. 4). SSA co-varies with AOD, similar 
to 6 March. However, the aerosol air masses in the study area 
are, in general brighter (higher SSA) on 15 March. Regions 
dominated by biomass burning aerosols (e.g., NW quadrant) 
might have different SSA values due to the burning of dif- 
ferent fuel types or the presence of biomass burning aerosols 
from different fire stages (e.g., flaming vs. smoldering) (Rei 
et al., 2005) on the two days. Aiken et al. (2010) mention that 
the fires were somewhat more intense in early March. Also, 
during the first 13 days of the field campaign (1-13 March 
2006) the surface conditions in MCMA were usually dry, 
with less than 10 % relative humidity (RH) in the afternoon. 
A cold surge during 14-15 March increased the RH to 40 % 
at the surface and 60 % at ^ 3 km (Fast et al., 2007) where 
smoke aerosols were found. Higher RH at the altitude with 
smoke (3 and 4 km) can also result in hygroscopic growth of 
the aerosols, contributing to higher SSA on 15 March. 

Based on the comparisons of MISR research retrievals 
against sub-orbital data, we confirm that MISR can distin- 
guish aerosol air masses qualitatively, based on shape, size 
and SSA. From MISR research retrievals, the following con- 
clusions are drawn about the two case study days. We see dif- 
ferent aerosol air masses on each day. On both days, spher- 
ical absorbing particles dominate over non-absorbing ones, 
although there is regional variation in this. The size of the 
dominant absorbing (SSA = 0.85) particles is smaller (r e ff = 
0.06 pm) on 6 March compared to 15 March (r e ff = 0.12 pm). 
Except for the central Mexico City area (SW quadrant), there 
is more dust on 6 March relative to 15 March. This result 
is confirmed by HSRL observations. Drier conditions on 6 
March favor the presence of more dust. Spherical absorbing 
particles having steeply sloping spectral absorption (smoke- 
like) are favored more than their flat (pollution-like) counter- 
parts on both the days. Smoke and pollution particles are dis- 


tinguished based on particle size. SSA from both AERONET 
and MISR research retrievals suggest that, on an average, the 
aerosol air masses are brighter (higher SSA) on 15 March. 

4.4 Analysis of MISR L2 standard aerosol product over 
our study region 

The MISR L2 standard AOD in the MCMA (our study area) 
was higher than the surrounding regions on both 6 and 15 
March 2006. AOD in MCMA is higher on 15 March than on 
6 March, all in agreement with sub-orbital and research re- 
trievals. However, the AOD on both days is lower than all the 
in-situ observations used in this study. AOD was not retrieved 
over a few regions within MCMA that resulted in either fill 
values (from neighboring retrievals) or no reported AOD val- 
ues. We further investigated the performance of the MISR 
standard aerosol retrieval algorithm over MCMA by look- 
ing closely at the retrieval diagnostic parameters and indica- 
tors reported for each 17.6 x 17.6 km retrieval region in the 
L2 Aerosol product. Successful aerosol retrievals identified 
primarily large and small/medium spherical non-absorbing 
particle mixtures. Mixtures containing dust were selected 
in few other instances. In one instance, the standard algo- 
rithm picked dust mixtures on 15 March in regions where 
thin cirrus clouds were present. When the retrieved AOD 
was very low (<0.05), a large number of successful mix- 
tures (>50) contributed towards the regional best estimate 
AOD. As discussed by Kahn et al. (2010), we find that in the 
urban environment of Mexico City, the standard MISR al- 
gorithm tends to pick mixtures with high single scattering 
albedo (SSA 558 = 1 . 0 ), resulting in an underestimation of 
AOD. The Regional Class Indicator flag shows that retrievals 
were not attempted over regions flagged as topographically 
complex. Retrieval flags (RetrAppMask) for regions with no 
L2 AOD retrievals also show that over most of these regions, 
retrievals failed due to poor correlation of the equivalent re- 
flectance spatial distribution from one MISR view angle to 
another. Poor angle-to-angle correlation can result from fea- 
tures such a complex terrain, aerosol plume edges, cloud 
edges, cloud shadows, and unscreened glint regions. We fur- 
ther investigated the diagnostic flags for six months of MISR 
L2 data for the same spacecraft paths (P026, P025) as the 6 
and 15 March observations, and find the above patterns to 
hold true for all L 2 retrievals over our study region having 
similar viewing geometry. The feasibility of retrieving AOD 
using the research retrieval algorithm and better agreement 
of AOD and aerosol properties with validation data indicates 
that the heterogeneity in aerosol particle type and aerosol 
loading over the 17.6x17.6 km retrieval regions of the MISR 
standard algorithm, the presence of very dark pollution parti- 
cles (SSA 558 = 0.7) in MCMA, complex topography and the 
vertical distribution of aerosols (smoke at ~ 4 km) limit the 
retrieval results from the MISR standard algorithm. This sug- 
gests the need to include darker pollution-like particles and 
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to perform higher spatial resolution retrievals with the MISR 
observations for urban regions such as Mexico City. 

5 Summary and conclusions 

In this study, we have assessed the ability of the MISR 
Research aerosol retrieval algorithm to distinguish different 
aerosol air mass types in the Mexico City Metropolitan Area 
on two days, 6 and 15 March 2006, during the INTEX-B 
field campaign for which we have coincident, ground-based 
and aircraft observations. Additionally, the MISR Research 
retrieval results were used to evaluate the performance of the 
MISR V22 Standard aerosol retrieval algorithm in our study 
area. Comparison of aerosol properties retrieved from MISR 
against the sub-orbital HSRL verify the ability of MISR to 
distinguish qualitatively between different aerosol air masses 
in MCMA. Both MISR and sub-orbital observations suggest 
the presence of mixtures of dust, biomass burning and urban 
aerosols in the study area. We summarize the findings from 
the MISR Research retrieval algorithm here. 

- In MCMA, MISR AOD shows regional variability at 
< 17.6 km scales on both 6 and 15 March 2006. The 
AOD and its variability are lower on 6 March than on 
1 5 March. Sub-orbital measurements confirm these re- 
sults and also provide qualitative verification of the 
MISR retrieved AOD pattern in the study area. Lo- 
cal meteorology governs the distribution of aerosols in 
MCMA. The basin-scale circulations on our case study 
days are found to be climatologically representative. 

- On the two study days, the MISR research algorithm 
cannot distinguish between weakly absorbing and non- 
absorbing particles. The weakly absorbing (<u 0 558 = 
0.95) spherical particles in mixtures were found to be 
replaceable by non-absorbing (&> 0 558 = 1.0) spherical 
particles in the retrievals. In regions where dust parti- 
cles were picked more frequently by the research algo- 
rithm, the dust particle (r e ff = 0.75 pm; <u 0 558 = 0.97) 
could be substituted by a larger non-spherical red-dust 
(ellipsoidal) particle (r e ff =1.18 pm) with a lower SSA 
(0.88). Based on earlier sensitivity studies (Kahn et al., 
2001; Kalashnikova and Kahn, 2006), this is expected, 
though this detail is specific to a complex urban envi- 
ronment such as MCMA. 

- In situations such as on 6 and 15 March 2006, when 
aerosols from local pollution mix with transported 
biomass burning aerosol from fires on adjacent moun- 
tains, the MISR research algorithm could not distin- 
guish between particles that show steep vs. flat spec- 
tral behavior of aerosol absorption or SSA (i.e., smoke 
vs. pollution). However, urban and smoke aerosols are 
distinguished based on particle size, as confirmed by 
suborbital validation data. 


- For HR retrievals in this study, obtaining robust re- 
sults from the research algorithm over the fire plumes 
to the east of Mexico City was difficult, because the 
plume appeared to move relative to the ground as the 
camera angle changed due to plume height, and the 
plumes were too small to co-register consistently at 
higher elevations. This also made the relative abun- 
dance of spherical vs. non-spherical particles diffi- 
cult to constrain at locations where discrete, elevated 
smoke plumes features were in the FOV of some of 
the MISR cameras. High-altitude cirrus clouds posed 
similar challenges on 1 5 March (as also discussed by 
Pierce et al., 2010). 

- Detailed analysis of MISR research retrieval results 
suggests that two distinct aerosol air mass types were 
identified based solely on particle microphysical prop- 
erties, over the Mexico City metropolitan area: one 
dominated by smoke and another consisting of a mix- 
ture of smoke, dust and pollution particles. In gen- 
eral, >70% of the AOD was attributed to spheri- 
cal absorbing particles and up to 15% to spherical 
non-absorbing particles. Non-spherical particle dust 
analogs contributed up to 30%. The gradients in the 
particle properties agree qualitatively with HSRL in 
most cases. 

In summary, over an urban area, surface heterogeneity, lack 
of specific knowledge about particle microphysical proper- 
ties, limited angular sampling, poorly constrained aerosol 
vertical distribution, coupled with low AOD conditions (as 
on 6 March) can make it difficult to constrain particle prop- 
erties. MISR research retrievals from this study demonstrate 
that over MCMA, MISR can distinguish at least two aerosol 
modes (fine and medium/coarse), and that fine mode (< 0.26 
micron) spherical absorbing aerosols are dominant in the re- 
gion during our study period. MISR retrievals cannot tell 
with confidence the difference between medium mode and 
coarse mode aerosol types in relatively low AOD situations 
over bright and varying urban surfaces. The distinction be- 
tween urban and biomass burning aerosols under these cir- 
cumstances is made based on particle size, though not on par- 
ticle SSA or its spectral slope, though such distinctions may 
be possible under more favorable retrieval conditions. Taken 
together, MISR was able to aptly distinguish regions dom- 
inated by dust, smoke and/or urban aerosols for our study 
days. Taking advantage of over 13yr of global observations 
from MISR, our future work will focus on mapping the sea- 
sonal and inter-annual history of aerosol air mass type and 
AOD gradients for mega-city environments globally. 


Supplementary material related to this article is 
available online at http://www.atmos-chem-phys.net/13/ 
9525/2013/acp-13-9525-2013-supplement.pdf. 
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